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A NEW NASTRAN CAPABILITY FOR DATA REDUCTION 

Michael Gallo and S . Mittal 
Bell Aerospace Textron 


SUMMARY 


A new mov* ~e, MODB, for the data reduction of NASTRAN results is described. 
NASTRAN analysis results can be filtered and sorted for minimum/maximum values and 
the printed output resulting from large NASTRAN runs can be limited based on a 
number of available user options, The sorting is done on stresses, forces and 
vector quantities like displacements, velocity, and acceleration. The module can 
be accessed via DMAP alters to existing rigid formats, and has been used on « 
large number of statics and dynamics problems at Bell Aerospace resulting in 
considerable savings in cost, time, and the amount nf printing. 


INTRODUCTION 


The high computational speed and large storage capacity of modern computers 
have enabled the analysis of large and complex structures. As a consequence, the 
structural engineer devotes much of his time to visually scanning, processing, and 
interpreting a larqe amount of finite element analysis results. This process is 
both time-consuming and error-prone. The way to alleviate the problem is to 
automate, whenever possible, the scanning and interpretation of the results, and to 
give the analyst the option to reduce the amount of computer output according to his 
engineering requirements. 

The need for these capabilities has been felt for some time and a number of 
proprietary post-processors have been developed by various NASTRAN users. However, 
none have been integrated into COSMIC/NASTRAN for general use. Bell Aerospace has, 
therefore, developed an engineering-oriei.ted module for data reduction which has 
been integrated into and offered to COSMIC/NASTRAN. MODB is currently in COSMIC/ 
NASTRAN release April, 1984, at Bell Aerospace Textron. No new Case Con ol or 
Bulk Data cards are required for input definition, only existing PARAM v..id DTI cards 
are used. 

This paper describes the data reduction module and its capabilities and demon- 
strates its application on an actual analysis. 


DESCRIPTION OF MODB 


MODB can process any real or complex OFP data block in SORTl format. This 
includes element stresses, element forces, loads, displacements, forces of SPC, 
eigem ^ctors, velocities, and accelerations. MODB currently contains coding for 
the 75 elements that existed in NASTRAN release April, 1982. The code is easily 
modifiable to include new OFP data blocks and elements as the need arises. 
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User requests are processed to update program controls used to suppress 
elements and to identify vector or element type components to be used for 
filtering and/or sorting. For complex vectors, stresses, or forces, the user can 
specify a real component mrntber, The vector component specified for either real 
or complex vectors must be in the range of 1 to 6. For both real and lex 
stresses and forces, the user input component number is checked to determine if 
it is in the range of valid component numbers for that particular element. The 
input data block undergoes a number of checks to insure it is valid for the current 
version of NOW. 

The stress, force, or vector quantities in the OFP data block are filtered 
and stored in compact form in memory for subsequent use. H* e packed data is 
sorted and efficient data reduction is performed in accordance with the u$er*s 
requirements. The filtered/sorted data block is output in standard OFP format 
with appropriate labels and printed by the existing OFP module. 

MODS may be run either in conjunction with NASTRAN as a single job or as a 
separate run following a regular NASTRAN analysis for which the necessary OFP 
data blocks to be filtered/sorted have been saved via CHECKPT or by using 
NASTRAN* s 0UTPUTT1 or 0UTPUTT2 modules. 


NOW CAPABILITIES 


The following is a brief description of the different capabilities presently 
available in HOW: 

A) Output Selection 

NASTRAN Case Control Set cards are used to define lists of node point lumbers, 
element numbers* or frequencies fer use in output requests in the normal manner. 
The NOW user can then limit the data reduction process to selective element 
types and components through Direct Table Input CDTI) if he so desires. 

B) Component Selection 

The user has the option to filter/sort element stresses and forces on any 
selected stress or force component. Vectors can also be filtered and/or 
sorted on a user specified component. In ail cases, a default component will 
automatically be selected by the program depending on the value of other 
fiiter/sort parameters input by the user. 

C) Sorting: Output in a Preferred Sequence 

Sorting on maximum magnitude, minimum magnitude, maximimi algebraic value, or 
minimum algebraic value is possible, as well as no sorting in cases where only 
filtering is desired. Multiple sorts in the same run may be requested. 

D) Filtering: Envelopes of Stresses, Forces, and Vector Quantities 

The user may define a value which defines a lower bound beyond which the 
search for sorting is to begin. Stresses and forces are filtered on a 



specified stress or force component* Vectors are initially filtered using all 
six components and then on the user selected component. All line items falling 
below the filter value are not printed. 

E) Reduced Output 

The user can control the number of Jims of printed output. A maximum of N 
lines will be output after filtering and/or sorting ha- been performed. 

MOOB input 


The main input to MOOB comes from NASTRAN OFP type data blocks generated by 
the regular analysis. The Case Control deck SET cards define the actual element 
identification ntmfeers and node numbers to be output for each subcase in the 
usual manner. The Bulk Data deck contains the main information required for 
data reduction by MOOB. This information is input through the use of existing 
PARAM and DTI cards unless defaults are used. 

In order to use MODS, the following DMAP statements must be included in the 
Executive Control deck either as a replacement for or an addition to the existing 
OFP module in any rigid format that uses SORT1 type OFP data blocks. 

MOOB OFPD, INDTI/OFPDX/C , Y , NUMCXJT = rf/C,Y,BIGER * /C , Y , SRTOPT»#/C , N , STRELYP*^ $ 

OFP OFPDX,,,, // S,N,CAFDNO $ 

The data block and parameter names used in the above DMAP statements rust be 
changed depending on the particular rigid format used and OFP data block to be 
processed by MOM. The Appendix contains a description of the input and output 
data blocks, and parameters used by MODB. 


SAMPLE APPLICATION 


A listing of the complete input and a sample of the output generated by MODB 
are given in figures 1 and 2. 


CONCLUSION 


Scanning of analysis results ami data reduction nay be performed with the 
aid of MODB in an automated fashion, thereby eliminating possible errors and 
waste of valuable man hours, both of which occur when performing data reduction 
in a manual and/or visual manner. As a result, the analyst is free to devote a 
larger portion of his time to engineering-oriented decision making based upon 
results obtained in an organized and comprehensive form. 

The use of MODB provides the analyst with an efficient and convenient tool for 
the study of KASTRAN analysis results and their presentation for project docu- 
mentation. 
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APPENDIX 


ENTRY POINT : MODB 

PURPOSE 

TO filter and/or sort S0RT1 QFP formatted data blocks. 

DMAP CALLING SEQUFNCE 

MODB OFPD, INDTI/OFPDX/C,Y ,NUMOUT=+0/C, Y,BIGER-0. 0/C, Y,SRT0PT=0 C,N,STRELTYP=0 $ 
INPUT DATA BLOCKS 

OFPD- — Any of the following OFP SORTl data blocks. 

a. Element forces (Major ID = 4 or 1004) 

b. Element stresses (Major ID = 5 or 1005) 

c. Displacements (Major ID = 1,1001,15 or 1015) 

d. Loads (Major ID = 2 or 1002) 

e. Force of SPC (Major ID = 3 or 1003) 

f. Eigenvectors (Manor ID = 7,1007,14 or 1014) 

g. Velocities (Major ID = 10,1010,16 or 1016) 

h. Accelerations (Major ID = 11,1011,17 or 1017) 

INDTI User input DTI which can control the elements and components to be 

sorted. INDTI may be purged. 

The DTI table consists of pairs of values B. and C. . 

ii 

B. - Element type identification number (Integer). If the B. value is 0, then the 

corresponding value is assumed to be a vector component. 

C. - Stress/Force component identification number on which sorting is to be per- 

formed (Integer). 

If the C. value is -1, the element type will be suppressed on the output file. 
An exampxe of this feature could be as follows: If an element type is to be 

sorted on two different values and output twice, this can be accomplished by 
two calls to MODB with two unique DTI tables* 

NOTES : 

a. Data block OFPD must be SOR'D, real or complex. 

b* If OFPD is purged or not recognized by MODB, then a non-fatal error will be 
generated and MODB will return. 

c. If INDTI is purged, the default sorting code will be determined by the value 
that satisfies the condition defined by parameter SRTOPT. 

d. INDTI can be used to modify the SORT codes as follows: 
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e. Each entry is on? of the following: 

1) For element types with existing data, two (2) words are used: 

a) Element type code (vectors use a zero). 

b) Stress item code. 

f„ The data item "ENDREC" must appear following the last word of the last entry 
input . 

g. A limit of 100 stress items may be handled. New element types having more 
than ICO items of stress data per element entry cannot be handled. The 
TRPLT1 elements have 65 items per element entry. 

OUTPUT DATA BLOCKS 

OFPDX Filtered and sorted OFP data block. 

Note ? OFPDX may not be purged. 

PARAMETERS 

NUMOUT lr.teger-input-default=0, NUMOUT controls the number of output lines. 

NUMOUT=0 implies that all items will be output after filtering and/or 
sorting has been done as controlled by parameters SRTOPT and BIGER. 

NUMCCJT=+N implies that only a maximum of N lines will be output after 
filtering and/or sorting has been done as controlled by parameters 
SRTOPT and BIGER. 

BIGER Real-input-default=0.0, BIGER is the filter value below which items 

will not be output. 

Since magnitudes are compared against BIGER, the default value of BIGER 
results in no filtering. 

Stresses and Forces are filtered versus a specific stress or force 
component. 

Vectors are filtered initially only if all six (real) degrees of 
freedom are less than BIGER and later, on only the component deter- 
mined by default or DTI input. 
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APPENDIX, con't. 


SRTOPT Integer input-def ault=0 . Controls the sorting option to be performed. 

Value Description 

-1 No sorting, 

0 Sort Oil maximum magnitude, 

1 Sort on minimum magnitude. 

2 Sort on maximum algebraic value. 

3 Sort on minimum algebraic value. 


STRELTYP Integer -input-def ault=0. Controls the element type to be processed 

for stresses and forces. 

Value Description 

C All element types will be processed. 

.GT.O Only element type STRELTYP will be processed. 
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PAUL 


NASTRAN EXCCU1IVF 


r o n i u u l 


DECK 


c h n 


ID DEMI 02, NASTRAN 

APP DISPLACEMENT 

SOL 1.0 

TIME 5 

ALTER 104 

HODB 0ES1 , /0ES1X/C , N. -1 /C ,N, 0.0/C, N,f)/C , N, 0 t 
OFP 0ES1X, , , , , // t 

MODE 0ES1 , /0ES2X/C, N, - 1 /C , N, 1 .0E2/C , N, 0/C , N, 0 $ 
OFP 0ES2X, , , , , // % 

ENDALTER 

CEMD 


SPHERICAL SHELL WITH PRESSURE LOADING. NO MOMENTS ON BOUNDARY FEBRUARY 11, 1983 RELEASE APR. 1984 

NASTRAN DEMONSTRATION PROBLEM NO. 1-2 

MIKE GALLO 

CASE CONTROL DECK ECHO 

CARD 

COUNT 

1 TITLE * SPHERICAL SHELL WITH PRESSURE LOADING, NO MOMENTS ON BOUNDARY 

2 SUBTITLE » NASTRAN DEMONSTRATION PROBLEM NO. 1-2 

3 LABEL - HIKE GALLO 

4 LOAD * 1 

3 SPC « 2 

6 OUTPUT 

? DISP « ALL 

8 STRESS - ALL 

9 BEGIN BULK 


»«* USER INFORMATION MESSAGE 207. BULK DATA NOT SORTED. XSORT WILL RE-ORDER DECK . 
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FIGURE 1 - EXECUTIVE & CASE CONTROL INPUT 



SPHERICAL SHELL WITH 1 
MS TRAN DEMONS THAT ION 

PRESSURE LOADING, 
PROBLEM NO. 1-2 

NO MOMENTS ON 

BOUNDARY 

FEBRUARY 

11* iy»5 RELEASE APR. 1984 

PAGE 7 

OES1X 

SORTED ON HA XI HUM MAGNITUDE VALUE , COMPONENT 7 , BIGER 

» 0.0 

, NUMOUT - -1 




S T ft E S 

S E S IN G 

ENERAL TRIANGULAR ELE 

(IN ELEMENT COORDINATE SYSTEM) 

M E N T S 

( C T R I A 2 ) 


ELEMENT 

ID. 

FIBRE 

DISTANCE 

STRESSES 
NOR HAL -X 

IN ELtMENT COORD SYSTEM 

NORMAL** Y SHEAR-XY 

PRINCIPAL STRESSES (ZERO SHEAR) 
ANGLE MAJOR MINOR 

MAX 

SHEAR 

26 

-1 .5O0QQOE+O6 
1 . 36000OE4O0 

3 . 1 41 ©28E402 
6.41 76486*00 

-1 .7587516*01 
2.7613416*00 

-1 . 6701346*01 
2. 370280E+01 

-2*8733 

42.9638 

3.1 4941 4E+02 
3.033722E+A1 

-1 . 8426356*01 
-2.11 7822E+61 

1 .6668396*02 
2.5747736*01 

24 

-1 .306666E+66 
1 .3000006400 

1 .6605726402 
9.2834326*01 

2.3373266*00 
-8. 445576E*01 

8.909364E*©1 
-4 . 4533046*01 

23.7115 

-13*3349 

2 .031 879E+02 
1 . 0339 17E+02 

-3.679332E*01 

-9.5613186*01 

1 .2099666402 
9.9262476401 

27 

-1 .30G000E+0© 
1 .366666E+66 

1 « 2796376*02 
-8.3888706*01 

4 . 6464756*01 
8. 175385E + 0© 

-5.3987876*01 
-1 . 61 3448E-01 

-26.4774 

-89.8993 

1 .3485456*02 
8.1 73639E*00 

1 .9573946*01 
-6.3868986*01 

6.7646276401 

4.6032326401 

28 

-i . soooooe+o© 
1 .3600601 >66 

1 .420986E*©2 
3. 413269E401 

1 .21 04946*01 
-1 .3552206*01 

-3.285495E*01 
2.6651 34E*01 

-13.4086 

23.5602 

1 . 4993066*02 
4.5741 1 26*01 

4.27301 0E 40# 
-2.714©44E*01 

7.2828786401 

3.6440896401 

23 

-1 .360606E+66 
1 .3000006400 

1 . 2531 43E*02 
-8.102733E*01 

-3 . ©371 49E*01 
1 .733096E*O1 

3 . 0005626*01 
-8.331900E401 

18.3460 

-60.2737 

1 . 443495E*02 
6 . 4902246*01 

-4.9606726*01 

-1.285906E402 

9.7078116401 

9.6730416401 

• 7 

-i.sooeeoE^oe 

1 .3000006400 

1 .2138776*02 
-8 .61 151 4E401 

-1 .323353E + G1 
5.02991 76*00 

-1 .4047476*01 
3.8204236*00 

-3*8102 

87*6552 

1 . 2281 70E*02 
5.1003546*0© 

-1 . 466290E*©1 
-0 .0271 56E*#1 

4.9739976401 

4.6728966401 

29 

-1 . 3000006400 
1 . 300000E400 

6.385297E401 
-4.921 7196*01 

5. 392976E*01 
-1 *21 72446*01 

-2.982939E+01 
1 .3698696*01 

-39*3490 
71 .8473 

9 . 031 0646*01 
-7.7103876*00 

2.947208E401 

-5.3479036*01 

3.0419286401 

2.2984226401 

16 

-1 .3000006400 

1 .3000006400 

6.3373126*01 

-7.5544336*01 

-2. 6079676*01 
1.585683E*01 

3.061 255£*01 
-8.271 609E+01 

23.9317 

-59.4603 

8.785620E401 

6.4637326*01 

-4.8502746*01 
-1 .2434306*02 

6.820947E401 

9.4501276401 

36 

-1.566606E+66 
1 . 5000006400 

0. 359991 F*01 
2.338O06E401 

1 .434454E*01 
-2.007843E*01 

-1 . 594507E*01 
8.715087E*0G 

-12*3624 

10.4190 

8 . 709440E*0t 
2.6994376*01 

1 . 084970E*01 
-2. 1 091946*01 

3.8122436401 

2.4343156401 

18 

-1 .366666E+06 
1 .5000006400 

6.91 43626*01 
-1 .0195856*02 

1.7413246*01 
-5 . 1606486*00 

-2. 131 494E+01 
6.021770E*©0 

-19.8716 

86.4338 

7.6939836*01 

-4.707470E*O© 

9.637009E*©# 
-1. 02331 66*02 

3.3651416401 

4.0772086401 

9 

-1.500000C400 

i .boeeeet+eo 

4 . 4648946*01 
-3.381 195E+01 

-4 . 2393786*01 
2 . 7388 756*01 

5.921 770E*G1 
-7.0786546*01 

26.8432 

-59.9183 

7. 461801E*©1 
0.0391036*01 

-7.2342B5E*01 

-9.4815©3E*01 

7.3490436401 

8.1603446401 

19 

-1 .5OOO00E400 
1 . 366O80E K>0 

6 . 1 ?6204E*01 
-3 . 359483E401 

1 .5010226*01 
-2.4©8O05E*01 

-2. 1 4091 3E*01 
2.391 499E*©1 

-21.2424 

60.8389 

7.00831 6£*01 
-1 .0735706*01 

6.68791 2E400 
-6.6939146*01 

3 . 1 49862E401 
2.8101736401 

4 

-1 . 500O00E400 
1 .5O0OO0E400 

4 . 1 380466*01 
*3 . 934A1 2E*©1 

-6 . 584657E+01 
4 .8032616*01 

5.5930946*01 
-5.221 292E*©1 

23.1059 

-64.9603 

4.524J88E*01 
7 .242371 E*01 

-8.970999E401 

-6.3737216*01 

7.747694E401 

6.6080466401 

15 

-1 . 500000E4©© 

1 .30O00OE40© 

3.577985E*01 

6.429749E-01 

-3.537155E+00 

-5.6?0733E*©1 

4.471141E*01 

-7.477538E+01 

33.1330 

-34.5094 

6 . 4963596*01 
3 . 2052846*01 

-3.2720926*01 
-1 .681 174E402 

4.B84225E401 

8.0085116401 

6 

-1 .5OOO00E4OO 
1 .3000006400 

3.484644E+01 

-3.348335E+01 

-4 . 6210336*01 
4 . 24491 0E + 01 

-3.478094E*01 
5. 7371 77E*©1 

-26.7525 
61 . 7475 

6 . 2461 29E+01 
7 . 327930E+01 

-7.38251 8E*01 
-6 , 431 355E*01 

6.8143236401 
6 .8796436 rOt 


FIGURE 2 - SAMPLE OUTPUT 
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0ES2X , SORTED ON MAXIMUM MAGNITUDE VALUE , COMPONENT 7 , B1GER * 100.000, NUMGUT - -1 

STRESSES IN GENERAL TRIANGULAR ELEMENTS < C T R I A 2 > 

(IN ELEMENT COORDINATE SYSTEM) 


ELEMENT 

FIBRE 

STRLSSES 

IN ELEMENT COORD 

• SYSTEM 

PRINCIPAL STRESSES (ZERO SMEAR) 

MAX 

ID. 

DISTANCE 

NORMAL -X 

NORMAL-Y 

SHE AR-XY 

ANGLE 

MAJOR 

MINOR 

SHEAR 

26 

-1 .5OO0O&E+0O 

3.I41028E*02 

~1 .758751E*01 

-1 .470134E+01 

-2.8753 

3.1 4941 4E+02 

-1 .842635E+01 

1.666839E+02 


1 . 50OO0OE+00 

0.41 7648E+00 

2 . 701 341E+00 

2.570280E+01 

42.9458 

3 .035722E+01 

-2.117822E+01 

2.576773E+01 

24 

**1 . 500000E+0O 

1 .66O572E+02 

2.33732AEt04 

8 9093A4E+01 

23.7115 

2.051 879E+02 

-3.679332Et01 

1 . 289904E+02 


1 .500000000 

9. 263432E+81 

-B.445574E*01 

-4.453304E+01 

-13.3369 

1 .03391 7E*02 

-9.501 318E+01 

9.920247E+01 

27 

-1 .500000E400 

1 . 279637E+02 

4. 646475E+01 

-5.398787E+01 

* 26.4774 

1 .548545E*02 

1 * 957394C+01 

6. 764027E+01 


1 . 5O06O0E+O8 

~8 . 388870E+O1 

8.1 75385E+00 

-1 .615448E- 01 

-89.8995 

8.1 75659E+00 

-6. 3BB898E+01 

4 .603232E+01 

28 

-1 . 5OOO00E+0O 

1 .420984E+02 

1 .2I0494E+01 

-3.2B5493E+01 

-13. 4080 

1 . 499304Et02 

4.27301 0C+00 

7.28287BE+01 


1 .500000E400 

3.41 5269E*01 

-1 .55522OE+01 

2.4651 34E+01 

23.5002 

4 . 5741 1 2E+01 

-2 . 71 4064E+01 

3.644089E«01 

25 

-1 .5000OOE+O0 

1 . 2531 43E+02 

-3.O37149E+01 

3.80Q562E+01 

18.3460 

1 .445495E+02 

-4.960672E401 

9 . 70781 1 E+01 


i .soooooEtoe 

-8. 102733E+01 

1 • 7330V6E+01 

-8.331 900E+01 

-60.2757 

6.490224E+01 

-1 . 285984E+02 

9. 675041E^01 

17 

>1 .5OOOOOE+00 

1 » 21 3877E+02 

-1 . 523353E+01 

-1 .404747E+01 

-5. 81 02 

1 . 2261 70E+02 

-1 .666290C+01 

6.973997E+01 


1.5O0O00E+00 

-8.81151 4E+01 

5.029917E+0© 

3. B20423E+00 

87.4552 

5 . 1 86356E+00 

-8.8271 54E+01 

4.472896E+01 
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S T R E S 

S E S IN G 

ENERAl. TRIANGULAR ELE 

(IN ELEMENT COORDINATE SYSTEM) 

M F. N T S 

( C T R i A 2 ) 


ELEMENT 

ID. 

FIBRE 

DISTANCE 

STRESSES 
NORMAL -X 

IN ELEMENT COORD SYSTEM 

NORMAL-Y SMEAR-XY 

PRINCIPAL STRESSES <ZERO SHEAR) 
ANGLE MAJOR MINOR 

MAX 

SHEAR 

1 

-I .500000E+OO 
1 .5QOO0OE+0Q 

-9.680287E+0© 
5. 1 441 41E+00 

-6 . 203754E+01 
4 . 7O4821E+00 

4 . 088985E^01 
-6. 393283E+01 

28.6858 

-44.9195 

1 . 26931 3E+01 
6.889757E+01 

-0.441 093E+01 
-5 .896B60E+01 

4. 855201E+61 
6.393309E+01 

2 

-1 .509O00E*©© 
1 .500O00E+00 

-9.91 2994E+00 
1 .066493E+01 

-3.837282E+01 

3.094138E+O1 

-6. 326425E+0 
3.6O9038E+01 

-38.6617 

52.8454 

4.070197E*et 
5 .829045E+81 

-0.B98776E+01 
-1 .66641 6E+01 

6. 484406E+01 
3 . 748730L+01 

3 

-i .5eooeoE+6e 
i .seoeooE+eo 

-1 .71554AE+01 
4 . 5902A3E+00 

-3 . 093681 E+01 
-4.491 492E+00 

5. 41 9623E+01 
-6.O51071E+01 

41 .3778 
-42.8542 

3.058640E+01 

6.073022E+01 

-7.867065E*01 
-6.0631 45E+01 

5.463252E+01 

A.OAseasE^ei 

4 
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